Fractal Dimension of Julia Set for Non-analytic Maps 
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z*^ + e are calculated perturbatively for small e. It is shown that Ruelle's formula for Hausdorff 
dimensions of analytic maps can not be generalized to non-analytic maps. 

I. INTRODUCTION 

The Julia set J of a map is the closure of the unstable periodic points 0-^] . It is an invariant set of the map and 
is usually a "repeller", that is, points close to J will be repelled away by successive iterations of the map. A simple 
example is the map on the complex plane: f{z) = z^, for which J is the unit circle. Points close to J will flow to one 
of the two stable fixed points: and oo. Thus J is the boundary or separator of basins of attraction. A much more 
complicated geometry appears for the Julia set of the map: 

f{z)^z' + c, (1) 

where c is a non-zero constant (see Fig. 1(a) for an example and Ref. Q for many other examples). In this case, J is 
a fractal and its topology undergoes drastic changes as c varies. 
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FIG. 1. The Julia set for (a) f{z) = z^ + e and (b) f{z) = z*^ + e, for e = 0.15 + i0.15. 

Before we proceed further, let us define a few notations. We denote /" to be n successive iterations of the map. 
That is /"(z) = f{f'^~^{z)). The set of all unstable cycles of length n is denoted by Fix/". Df is the derivative 
matrix of /. If / is an analytic map, i.e. du/dx = dv/dy and dv/dx = —du/dy with /(z = x + iy) ~ u -\- iv, then 
detDJ ^ \df/dz\^. 

For analytic maps, the Hausdorff dimension Dh of the Julia set J can be calculated with a formula due to a theorem 
of Ruelle ||: 

lim A„(Dff) = 1, (2) 

n — '■oo 

where 

Using the formula, Ruelle and Widom et al. calculated Dh for the map (|l|) in powers of c for small |c|. It was 
not clear then whether the formulas (|^) and (|^) can be generalized to non-analytic maps. The natural generalization 
of (0) to non-analytic maps would be 



1 



A^{D)^ J2 \det or \-^/\ (4) 

The calculations I present below show that the combination of ^ and (^) does not give the correct Dh for non-analytic 
maps in general and Dh can be calculated directly with the perturbation theory developed in Ref. [p[. 



II. THE MAP f{z) ^ + ez* 
Let us first consider the non-analytic map 

f{z) = z' + ez*, (5) 

where * denotes the complex conjugate. When e = the Julia set J is the unit circle and can be parametrized as 
z{t) = e^'^**. The map on J is 

f{z{t)) = zi2t). (6) 

When e ^ but small enough so that J is topologically equivalent to a circle we can still parametrize J so that Eq. (||) 
is satisfied |^,^ . If a map /c with a parameter e satisfies 

[Mz)Y^Uiz*), (7) 

then 

zeFixf^ ^ z* eFixf^,, (8) 

which implies that 

J(/,) = [J(/,.)]*, (9) 
where J(/) is the Julia set of /. In particular, if J can be parametrized as z{t) then 

z,{t) = z:4~t). (10) 

It is easy to see that the map (||) satisfies Eq. (^. 

Following Widom et al. ||^, we formally expand z{t) in powers of e 

z{t) = e2"*[l + eUi{t) + e*Ui{t) + e^U2{t) + e*^U2{t) + ee*U2{t) + ■■■], (11) 

where the functions Ui{t),Ui{t),U2{t),U2{t),U2(t), . . . are all periodic with period 1. Eq. ( p^ ) implies that all the 

functions U{t) satisfies U{t) = U*{—t). Substituting (|ll| ) into (^ and equating terms with the same power of e, we 
get 

C/i(2i)-2[/i(t) = e-6"*, (12) 

f7i(2i) - 2C7i(<) = 0, (13) 

U2{2t) - 2U2{t) = Ul{t) + e-S"*f7i*(i), (14) 

1/2(20 -2[72(t) = c7i'(t), (15) 

Wt) - 2U2{t) = 2Ui{t)Ui{t) + e-^""f/*(t). (16) 



The solutions are 
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= 0, 




oo - 







g-37ri2^(2'="i+2'"i)t 



(18) 

2j+k+l ' ^^^^ 

(20) 

^ITi ■ (21) 
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It is easy to see from Eq. (||) that unstable cycles of length n are 

Fix/" = {z{t,) : t, = = 0, 1, . . . , 2" - 2}. (22) 

We now evaluate A„{D) as defined in (|^). Note that 

det D r = if det ( f ^ + f 'J'! 

11 \2j/i + Im(e) 2a;i - Re(e) 

I 12 »-l 



m=0 

where the last equality holds to the second order in e. Denote 

2"-2 



where tj 's are given by Eq. (p^. 

E I det Drr^/^ 



< >n= 2^7^ E ^(^j)' (23) 
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Substituting Eqs. (hj)-(ml) into (O) and using the identity 



^ 2TTimt ^ _\ 1, m = 0mod2"-l , . 

^ 10, m 7^ mod 2" - 1 ^^^^ 



it can be shown, after some algebra, that 



|z(2'%.)r^>„=l + |ep(^ ^_ + ), (^>2). (26) 



m=0 

Substituting (|i) into (|2|) yields 



A„(i?) = 2"(i-^)[l + |.p(^-^)], (n»l). (27) 



If we were to use Eqs. (27) and (Q) to obtain a Hausdorff dimension, we would get Dh = 1 — kl /(81n2), a value 
smaller than 1 for small but nonzero e. We show in the following that this value of Dh is incorrect. 
Let 

where z{tj) 6 Fix/" (Eq. (|2|)). The Hausdorff dimension Dh of the set Fix/" in the limit n ^ oo is such that 

lim XniDn) = 1. (29) 

n — >oo 

This Dh should also be the Dh of the Julia set J. We now evaluate Xn{D) to the second order in e. Putting 
Eqs. (0)-(|l|) into Eq. (|l|), we write 
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z(t,+i)-z(t,) = Co + Ci|e| + C2|er 

Then to the second order in e, 

XniD) = lgp(2" - 1){1 + 1^ < Re(C*Ci) >„ 

+ ^ < ICiP >„ + < Re(Co*C2) >„ +^ < (Re(Co*Ci))2 >„ 
where Eq. ( p3|) is used. With the help of the identity (^) we get 

27r 

ICol^ = 2(l-cos^^), 

< Rc(Co*Ci) >„ =0, (n > 2) 
< >„ =F(n), 

<Re(Co*C2)>„ =^(1 + ^), 

< (Re(Co*Ci))2 >„ = i|CoP < >„, 

where 

N 2 1 „ 3 • 2'=-! - 1 

Fn = --2> ^cos27r . 

^ ^ 3 4fe 2" - 1 

fc=i 

The function F{n) can easily be solved for large n in the following way. Note that for n >> 1 

^, s 1 / 3 • 2t7 , 1 „ , s 

^^('^ + 1) = 2 ~ 2^) + 4^(") 

= — rr + -i^(ri). 

4«+i 4 ^ 

Substituting F(n) = i/(n)/4" into Eq. (|8|), we have 

Hin + 1) ^9tt^ + H{n), 

which has the solution 

H{n) = 97r^n + a, 



where a is a constant independent of n. From Eqs. (|3l|), ( |32[ ) - (36), and ([40[), 

16 

Eqs. @ and imply 

1 + — ^|e|2. 
161n2' ' 



III. THE MAP f{z) = z*'^ + e 



Next, we consider the non-analytic map 



The map (Hj) has the property of Eq. (0), so that J{f^) = [J(/e-)]*- 
Let us parametrize J in such a way so that 

fiz{t)) = z{-2t), z{t) e J. 

The set of unstable cycles of length n is 

Fixr = {z(t,):i, 



(44) 



j = 0,±l,±2,...}. 



(-2)" - 1 

The number of elements in Fix/" is |(— 2)" — 1|. Following similar procedures as in the previous section, we have 



(45) 



Ui{t) 
Ui{t) 
U2{t) 

U2(t) 



oo 

fe=l 



4fc 



4fc 



OO 
OO 



4i+fc+i 



g~-7ri4^(4'""-'+4'"-')t °o g-7ri(4'=+4')t 

-+ E 



-7ri2^(4'=+4V2)t 



fc,i=l 



4fc+i 



2J+2k+2l 



AniD) (Eq. (|)) and Xn(£') ((Eq.(|28l)) can be calculated to be 



Ar^iD) = Xn{D) = 2"(i-^)(l + \nDW\ 



In this case, An{D) ~ Xn{D) and it gives the correct Hausdorff dimension 

l^|2 

Dh = 1 
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(46) 
(47) 
(48) 
(49) 
(50) 

(51) 
(52) 



The reason for Ruelle's formula to work in this case is that for the non-analytic map (^3|) / (z) is analytic: 

P{z) = {z^^e*f + e, (53) 

and that J{P) — J{f)- Note that ( |5^ is the same as the Dh of the analytic map (|l|) f{z) — z"^ + e |^,^, to the 
second order in e. Indeed, /'^(z) and thus J are identical for the two maps (^ and ( ^ ) for real e. For complex e, 
however, the two Julia sets look quite different (Fig. 1). 



IV. DISCUSSION 

Since Ruelle's formula relies on the analyticity of the map, it is no surprise that it brakes down for non-analytic 
maps. When J is a closed curve, Dh can be calculated from Xn{D) (Eq. (^)) for both analytic and non-analytic 
maps. When J is no longer topologically a circle, it can be difficult to utilize a formula based on distances between 
unstable cycle elements. In this case, it remains a challenge to formulate an efficient method for the calculation of 
Dh for non-analytic maps. Finally, the quantity An{D) (Eq. (|^)) can be very useful even for non-analytic maps. For 
example, it can be used to calculate the escape rate for points close to J [Q. 
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